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A general type of mathematical argument is described, which applies to all the cases in which 
dynamo maintenance of a steady magnetic field by motion in a uniform density is known to be im- 
possible. Previous work has demonstrated that magnetic field decay is unavoidable under conditions 
. of axisymmetry and in spherical or planar incompressible flows. These known results are encom- 

' passed by a calculation for flows described in terms of a generalized poloidal-toroidal representation 

of the magnetic field with respect to an arbitrary two dimensional surface. We show that when the 
velocity field is two dimensional, the dynamo growth, if any, that results, is linear in one of the 
f"^ . projections of the field while the other projections remain constant. We also obtain criteria for the 

existence of and classification into two and three dimensional velocity results which are satisfied by 
^ ^, a restricted set of geometries. In addition, we discuss the forms of spatial variation of the density 

I ' and the resistivity that are allowed so that field decay still occurs for this set of geometries. 
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C/3 i 1. Introduction 

a 

The steady state dynamo problem in a uniform density can be stated as follows. Given a uniform 
" electrically conducting fluid, contained in a volume, what condition of must be satisfied by the 

. velocity v of a steady motion in order that a steady field B can be maintained by dynamo interaction 

' between the motion and the field? This classical problem of the magnetic dynamo concerns the 

question of the amplification or maintenance of the magnetic field in cases where the induction 
equation is valid. The equations yield only decaying solutions when the velocity and magnetic fields 
are both axisymmetric [I, 2] or if the geometry has planar symmetry [3]. In both situations, the 
velocity and magnetic field can be three dimensional but do not depend on at least one of the 
coordinates. In situations where the velocity is two dimensional, but the magnetic field is three 
dimensional, the impossibility of dynamo action has been proven if the flow is planar [4, 5, 6] or 
spherical [7, 8]. 

In the following, we use a generalized toroidal-poloidal representation of the magnetic field, B, 
with respect to an arbitrary two dimensional surface and derive two scalar equations for the poloidal 
and the toroidal potentials from the induction equation. We also prove a result, which is an extension 
of the one in [9], that incompressible two-dimensional velocity flows in situations other than in the 
above mentioned antidynamo theorems (where the field decays) , lead to linear growth in one of the 
field components and are otherwise slow. Here, the fluid velocity is two dimensional in the sense 
that it lies entirely on surfaces which can be described by x(r) = constant. The approach taken 
here also lends itself to a unified and simpler exposition of the previously cited results. In addition, 
we consider the special cases of spatially variable forms of the density and magnetic diffusivity. 
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A closer study of the cases in which the dynamo maintenance of a steady field is impossible 
may throw light on the general dynamo problem. The results suggest that the number if cases in 
which dynamo maintenance is impossible is restricted. We derive such criteria for the existence of 
antidynamo result for a given geometry. 

2. Normal projection of the induction equation 

The starting point of all the above investigations is the well known induction equation in MHD 

QtB + V X (B X v) = -V X (?7V X B), (1) 

with the constraint 

V • B = 0, (2) 

where r] = c^/Ana is the magnetic viscosity and a is the conductivity. As in the above cases, the 
fluid is assumed to be incompressible ( V • v = 0) . We make the additional simplification of taking 
rj to be uniform. Later, we discuss the effects of relaxing the latter simplification. 

The RHS of the equation (1) represents resistive dissipation whereas the LHS contains a term 
which represents stirring of the B field by fluid motions. In analogy with the heat conduction 
equation, this shows that in a static fluid the fields decay, while stirring may induce field generation. 

Consider the fluid velocity and the magnetic field to be described in terms of components 
perpendicular and parallel to the surfaces defined by x(r) = constant. The quantities, = B • Vx 
and Vj(. = V • VX; satisfy the reduced form of the induction equation projected normal to the surface 

dtB^ - B ■ Vv^ - TjiV^B)^ = 0, (3) 

where dt = dt+v-V, and the constraint (2) and incompressibility condition was taken into account. 
One can expand the third term of (3) using identities (35)-(37), and obtain 

dtB^ = r?(V2Bx - V • [(Vx • V)B + (B • V)Vx]) + B • Vv^, (4) 

where the last term, B • Vv^, equals zero when the velocity fields lie on the surfaces, x(r) = 
const. Now employing the identity (36) on the LHS, writing the first term on the RHS in terms of 
V ■ (B-^VB-^) after multiplying throughout by B^, and integrating, while taking v^ = (we relax 
this later), this can be further reduced to the dissipation theorem [9], 

idt ^ B^d^r = -r?^{(VB^)2 + 6^9(3, x)}d^r, (5) 

where a surface integral over \7 -{B^VB^) obtained from Gauss's theorem vanishes at large distances. 
Here we have introduced a useful quantity 

e(Y,x) = V • [(Vx • V)Y + (Y • V)Vx] 

= [25fe(y,)+>fea,]9Ax + Vx- V(V-Y). (6) 

It can be seen that the last term in equation (6) vanishes if Y is solenoidal or if V • Y is independent 
of a coordinate directed along Vx- The only positive contribution, leading to growth, can come 
from the tensor term in O on the RHS of (6). Clearly, if the surface is planar (x(r) = z) or spher- 
ical (x(r) = J'^/2) the term becomes zero {didk equals or Sik, respectively, and in the latter case 
the condition, (2) needs to be further applied), implying that B-^ decays. This point was made by 
[9, 10]. Hereafter, we suppress the notation 6(Y,x) to 6(Y) unless x is specified. In the following. 
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where we keep the treatment general (by keeping the velocity and magnetic fields three dimensional), 
we find the same term occuring in the surface projection of the induction equation. This enables 
us to expand upon these conditions for field decay and generalize the antidynamo results cited earlier. 

3. Parallel projection of the induction equation 

It is convenient to express the magnetic field in terms of the local "poloidal" and "toroidal" compo- 
nents 

B = Bp + Bt = VxVx (V'Vx) + V X ($Vx), (7) 

where il> and $ are the generalized poloidal and toroidal flux functions. This is analogous to the 
description of magnetic field given by [11, 4] for spherical geometry. We express the field in a 
coordinate system in which a special coordinate, q, is normal to the surfaces and is given by 

X(r)=/(g). (8) 

The connection to the corresponding formulae in spherical geometry lies in the fact that any smooth 

surface has a local radius of curvature. Therefore, wc implicity demand that the surface be smooth 
(have first derivatives defined). Next, using properties (39, 34) it can be seen that 

= V X (V^ X Vx) • Vx 

= V • [{Vi) X Vx) X Vx] 

= V • [(W • Vx)Vx - (Vx)'VV'] 

= -(Vx X V)V = -VfV. (9) 

In order to examine the local components parallel to the surface, we take a normal projection of the 
curl of the induction equation (1), 

dt[{W X B) • Vx)] + Vx • [V X V X (B X v)] = -Vx • [V x V x (r/V x B)], (10) 

After expanding B, the operand of the time derivative in the resulting equation can be reduced 
in a fashion similar to Eq. (9) 

(VxB).Vx = -Vf$ + e(Vx [7/.Vx]). (11) 
Defining, C = V x B, and using Eqs. (35) and (11), the term on the RHS of Eq. (10) yields, 

-7?(curl^B) • Vx = r/(V^C) • Vx 

= r, (V^Cx - V ■ [(Vx • V)C + (C • V)Vx]) 

= -77V2Vjj$ - r/e(C) +r/V2e(V X (7/;Vx)), (12) 

after performing manipulations identical to those required in obtaining Eq. (6). The second term 
of equation (10) after applying Eq. (39), is 

-Vx ■ {V X V X [v X (V X V X (^/-Vx)) + v x (V$ x Vx)]}. (13) 

The second term may be evaluated in steps as follows 

V X (V$ X Vx) = v^V* - (v • V$)Vx 
V X [v X (V$ X Vx)] = Vvx X V$ - V(v ■ V$) x Vx 
-Vx ■ V X {V X [v X (V$ X Vx)]} = -V-[V$(Vx-Vv;,)- Vv;,(Vx-V$)] 

-V^(v-V$). (14) 
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Defining, P = v x Bp, via (A.l) reduces to 

I?x = V (Bp X Vx). (15) 
The first term in Eq. (13) will reduce to 

-Vx • [V X V X D] = -e(X)) + V^V^ (16) 
where the properties in the calculation (12) are used. 

4. Results 

Now one can write the cqiiation describing the evolution of the parallel components of the field 
by including all the terms simplified to the forms given in Eqs. (9) (15), and rearranging terms, as 

V|(dt$)-r/V'v|$-r/e(C) = V • [Vv;,(Vx • V$) - V$(Vx • VvJ] 

+ {dt - r?V2)e(V X [i^Vx]) - @{V) + V^V^. (17) 
The full form of the corresponding equation for the normal component is 

[dt - r/V2]V^V - '?e(B) = -(B . V)v;„ (18) 
where the second term involves the normal component of velocity. 

4-1 Non-diffusive flows (rj = 0) 

An exclusion theorem proved by [9] that states that two dimensional non-diffusive flows with a 
stationary velocity field and the property v • a = (the flux helicity density) under the gauge 
condition, V • a = 0, where a is the vector potential for the velocity, will lead to a conservation of 
B • a. We sketch their proof below. We write the Euler equation as 

^tv = f ; f = Vw - (v • V)v, (19) 

where all potential forces are collected in w. We can then write, under a gauge condition V ■ a = 0, 

dt(a • v) = V • V(a • v) + 2a • f + V • (a X curP^f). (20) 

Next we substitute for f and use 

-2a((v • V)v) = V • (au^) + 2(v x w) • a (21) 

enroute to obtain 

dt(v ■ a) = 2a • (v X w) + V • [(v • a)v + a X curP^f + 2(w - v2/2)a]. (22) 

This implies that the flux helicity of the flow lines, Hy = J v • a d'^.x, is conserved for Beltrami flows 
(v (X Lo), or for potential flows {lo = 0). Similarly, in the flux freezing limit of the induction equation, 

QtB = V X (v X B); dta = curP^f (23) 

we can obtain the following after some transformations 

dt(a • B) = B ■ V(v • a) + B • curP^f. (24) 
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Now it is easy to see that for stationary flows (f = 0), if the flux hehcity, Hy, is zero everywhere, 
that the cross helicity, i? c = / B • a d^x, is steady. 

We now generalize this result by not imposing any limitation on the gauge transformation- as it is 
impossible to always simultaneously satisfy the conditions v • a = and V • a = 0. When the velocity 
is two dimensional (v^ = 0) and r/ — > 0, it is clear from (18) that is an integral of motion and can- 
not grow with time. Further, if one takes a stationary velocity field with v ■ a = and V • a = 0, this 
implies that the flow lines are along the intersection of two surfaces; i.e. v = x Vx- As a result, 
is another constant of motion (v| =0). So the induction, v x B = B^V^ — B-^^Vx, is independent 
of time and the growth of the remaining component of B, (along v), can utmost be linear. Hence, 
we can conclude the intuitive result that flows with zero linkage (and V • a = 0) cannot lead to an 
exponential growth of the field, but utmost to a linear growth of the field in the direction of the flow. 

4-2 Antidynamo theorems 

Based on the equations (18) and (17), we can immediately divide the antidynamo results for two 
and three dimensional velocity fields. For the 3D results, we take both the velocity and the magnetic 
field to be three dimensional but invariant along the special coordinate q {dq = 0). In the case of 
2D results, the velocity field is two dimensional, v^ = 0, and the magnetic field is three dimensional. 
Under either of these conditions the first term on the RHS of (17) vanishes. The advection-diffusion 
operator (dt — r/V^) can only manipulate the field and cannot cause growth (cf. (5)). 

Here, an antidynamo case is defined as a situation in which the flux functions, tp and <I> ^ 
everywhere (as t oo) with the boundary conditions that these flux functions vanish at remote 
surfaces which enclose the volume of fluid. So the strategy in finding antidynamo situations is to 
identify the conditions when one equation completely decouples from the other and the decay of the 
corresponding flux function kills the source terms in the other. It is natural to consider spherical 
and planar geometries first, since 8(B) and 0(C) would then be zero (cf. (6)). It is to be noted 
that 0(V X [^Vx]) is zero for spherical, planar and axisymmetry. 

4-3 The cases of spherical and planar geometries 

For spherical geometry (x — ^r^) or planar geometry (x = z), the equations (17)-(18) reduce to 

[dt-r/V2]$ = Q(V,v) (25) 
[dt - T?V^]VjjV = -B • Vvx (26) 

where 

Q(^, v) ^ [Vlr\v^V^ - 9(1?) - V • [V$(Vx • Vv^)]}, (27) 

and V| represents , the angular momentum operator in the spherical case or + dy in the planar 
case. In the two dimensional case (v^ = 0), we have a source term only in the toroidal equation, (25), 
while tp decays. According to the definition of 2?. as V' ^ 0, P ^ and the RHS of (27) vanishes. 
Therefore Q ^ and $ will decay when t ^ oo. This follows from the arguments after (5), and from 
potential theory which demands that the mean value of [V|] (0) is zero in the volume enclosed 
by a surface on which it vanishes. Physically, the normal component difi^uses out and the fleld is 
confined to two dimensions; as a result, the field is transported like a scalar, and hence decays. 

Now, in the three dimensional planar case (x = -2,52 = 0) there is a source term only in the 
poloidal equation, (26). It can be easily seen from (6) and (15), that the toroidal source terms 
involving T> and Vx • V = are zero. Subsequently, as $ — > 0, B • V — > Bzdz{= 0) and ip decays. 
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It is interesting that in the three dimensional spherical case none of the source terms in the above 
pair of equations are zero to begin with, and hence dynamo action occurs. 

4-4 The case of axisymmetry 

In axisymmetry {x = 07 • V = w~'^d^ = 0) and V is poloidal. As a result, Q = 0, as the 
term involving V is zero [cf. (16)], while 0(A, </)) = —{2/w)dj^{A^/w), expressed in the cylindrical 
coordinates, {zu,(f), z). Then the equations, (17)-(18), simplify to 

Z)2[dt-,yD2]$ = 0, (28) 

[dt - T]^-^^^w^]{B^/vu) = -B • V(v^/tu), (29) 

where = — -^d^u, and is known as the Stokes operator, and V| represents [Ti7~'^(j) x V]^ = 
zu~'^D^. Since the term can be reduced to a divergence term and subsequently to a surface 
integral by Gauss's theorem which vanishes due to the dipole behavior of the field at large radius, it 
does not contribute to growth; see for example, [4] (p. 114). Therefore $ ^ 0, B • V — >■ B^d^{= 0) 
and tp decays according to similar arguments in [4] ( p. 115). 

4-5 Effects of spatial variation of density and resistivity 

It is easy to sec from the form of the continuity equation in steady state or under the anelastic 
approximation when the density is dependent only on the special coordinate q, 

pV • V + VgVgP = 0, (30) 

that V • V = is valid when v, = 0. Therefore only the above 2D results still hold. Non-uniform 

resistivity introduces a term (Vrj x Vx) ■ (V x B) in the equation (18) for ip. This term is zero if 77 
is a function only of q and hence does not alter the decay of ip. This was commented upon by [6] 
for the cases of spherical and planar geometries. We now show that $ also decays in some special 
cases. The spatial variation of 77 in q, however introduces a nonvanishing term {dzi])(dz^) for planar 
geometry or {l/r){drr]){dr[^r]) for spherical geometry in equation (25). The 3D planar case follows 
trivially. Now when v^^ = 0, one can invoke a theorem on the resulting elliptic equation [12, 5], 
which states that only the constant solution ($ = 0) is possible under the condition of $ vanishing 
at large distances. Similarly for axisymmetry, r]{4>) introduces the term {d^r])C^ in (28). When 
there is no differential rotation, — > and which depends solely on B^, vanishes and $ decays 
as before. Therefore, the above 3D planar and 2D planar, spherical, and axisymmctric results arc 
still valid if 77 is a function only of q. Also, if there is no differential rotation, axisymmetric fields 
cannot be maintained if ri depends on </>. It shown that in axisymmetry [13], poloidal fields cannot 
be maintained even by a compressible fluid with 77 as a function of space and time. 

5. Concluding remarks 

In this paper we have cast the induction equation (17)-(18) in a geometry given by the surfaces 
X(r) = constant. This was useful in extending a previous result for incompressible two dimensional 
flows while unifying, classifying, and simplifying the proofs of the known antidyanamo results, and 
thereby providing some new insights into the structure of the induction equation. In order to deduce 
the general conditions of decay, we can consider the order of decay of the flux functions, and 
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If $ is to decay first, then all the potential source terms which are on the RHS of (17) containing tp 
and V should be zero. This includes the condition that 



(31) 



The above equation automatically ensures that 0(V x [V'Vx]) = Vx • [V x V x V x (V'Vx)] is zero 
(c.f . (16)). The equation, (31) is true for planar and axisymmctry. Consequently, the 3D results of 
planar and axisymmetry follow, as the first term on the RHS of (17) which involves Vx • V is zero. 
On the other hand, if ip is to decay first then the RHS of (18) should be zero demanding the 2D 
condition, y-^ = 0. Further, 



must hold to ensure that tjj decays as per (5) and the boundary conditions that the flux functions 
vanish at the remote boundaries. This takes care of the decay of $ since the remaining source terms, 
which involve tp, in (17) vanish. The condition, (32), is true for planar, spherical and axisymmetric 
geometries. The unifying aspect of treatment of the boundary conditions used here is that [Vjj]~^(0) 
is zero with vanishing flux at the boundaries. 

We have been able show the impossibility of dynamo action (as defined in §4) with rj ~ const 
and V • V = 0, exists only for a restricted group of geometries as allowed by (31) and (32) for the 
case of 3D and 2D results respectively. In addition, the above results need further qualifications that 
V is bounded for all time and that "spiky" time dependent behavior is excluded [14]. 

As suggested in a classic paper [15], using a different approach, that the impossibility of a 
steady dynamo for a given geometry depends on the existence of an arbitrary current, j', such that 
/ j • j'd'^r = 0, and that only a restricted set satisfies this equation. Here we have derived specific 
conditions that determine these geometries. However, a more rigorous analysis of the above two 
conditions is needed to find the set of all possible x that satisfies the above criteria. In a paper, in 
preparation, I investigate the geometries that satisfy the above criteria. 

Acknowledgment: I thank P. H. Roberts for his valuable comments and for a critical reading of 

the manuscript. 

Appendix: Formulae referenced in the text 




(32) 



a • (b X c) 
a X (b X c) 

V X (V X a) 



b • (c X a) = c • (a X b) 
b(c • a) — c(a • b) 

V(V ■ a) - V^a 

(a • V)b + (b • V)a + a x (V x b) + b x (V x a) 
b • (V X a) - a • (V X b) 
a(V • b) - b(V • a) + (b • V)a - (a • V)b 
VV'{x,-}a + VV{x,-}a 



(33) 

(34) 
(35) 
(36) 
(37) 
(38) 
(39) 



V(a-b) 



V • (a X b) 
V X (a X b) 
V{x,.}(^a) 
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